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Helioseismology is the study of the variations in the internal structure and properties of the dynamics of the Sun from mea- 
surements of its surface oscillations. With the 2010 launch of the Solar Dynamics Observatory (SDO) we are undoubtedly 
approaching a new dawn for local helioseismology, as the extent and quality of raw surface oscillation data has never been 
better. However, advances in theory and modelling are still required to fully utilise these data, especially in magnetic active 
regions and sunspots, where the physics is poorly understood. 
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1 Introduction 

Understanding the origin of the Sun's magnetic field is ar- 
guably the most important topic in solar/stellar physics to- 
day. Sunspots and starspots are the most visible manifesta- 
tions of solar/stellar magnetic activity, with sunspots tend- 
ing to appear at well-defined latitudes on the Sun, varying 
with the 11 -year solar cycle. Our comprehension of, and 
capability to predict, the 11 -year activity cycle of the Sun 
is far from being complete however; as suggested by the 
recent unusually long solar minimum of Cycle 23. Under- 
standing the mechanism of the solar dynamo is therefore of 
the utmost importance, not only to be able to explain this 
collective behaviour in the Sun, but to also explain the fun- 
damental physical processes of the space environment from 
the Sun to Earth, to other planets, and beyond to the inter- 
stellar medium. 

Sunspots can provide us with important clues on how 
the solar dynamo operates, therefore a comprehensive dis- 
cernment of their subsurface evolution, structure and dy- 
namics is essential in order to establish accurate physical 
relationships between internal solar properties and magnetic 
activity in the photosphere. Helioseismology is the only tool 
available that can probe the structure and dynamics of the 
Sun through the study of acoustic oscillations. 



1.1 Local Helioseismology 

Helioseismology is a diagnostic tool that allows us to probe 
the solar parameters to compare with theory and obser- 
vation. Over the last fifteen years, considerable progress 
has been made in the understanding of the solar interior, 
mainly as a result of the precise determination of the fre- 
quencies of millions of individual global modes of oscilla- 
tion using observations from the Michelson Doppler Imager 



(MDI; IScherrer et al.|1995| ), on board the Solar and Helio- 
spheric Observatory (SoHO) spacecraft, and ground-based 
networks of telescopes. Measurements of flows and temper- 
ature variations in the solar interior have already led to an 
increased understanding of solar activity; for example, the 



detection of subsurface meridional flow (e.g., Giles et al. 
,1997^) has had a major impact on theories of the solar dy- 
namo. 

Local helioseismology has also shown that wave prop- 
agation is different in sunspots than in the quiet Sun ( |Gi-| 
|zon et al.||2009t [Moradi et aL]|2010| l, with the wave travel- 
time shifts in the vicinity of sunspots typically interpreted 
as arising predominantly from magnetic fields, flows, and 
local changes in sound speed. The magnitude of the travel 
time shifts also appear to acutely depend on the details of 
the data analysis, the wave travel distance, and frequency 



(|Braun & Birch|2008t |Couvidat & Raj aguru|2007| [Moradi 



et al.|2009[ [Thompson & Zharkov|2008[ )T 



Various local helioseismic techniques exist to analyse 
the oscillations. For example, Fourier-Hankel analysis has 
been used successfully in the past for studying the interac- 
tion of acoustic waves with sunspots (Braun et al. 1987[ l. 
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This process decomposes the solar oscillation signal, ob- 
served in an annulus around a sunspot, into inward and out- 
ward propagating wave modes. Another technique is the 
ring-diagram method ( |Hill|[1988 ). Ring-diagrams are local 
spatial power spectra at constant frequency. Cuts at constant 
frequency through the 3D power-spectra reveal nested el- 
lipses or rings that change shape and shift centre under the 
influence of alterations to the solar medium. Time-distance 
analysis is another method which computes the travel time 
of a wave packet travelling between two points on the so- 
lar surface ( [Duvall et al.|1997 1. Finally, there is the method 
of helioseismic holography (acoustic imaging), which re- 
constructs the subsurface wave fields by propagating waves 
either forward or backward in time ( [Chang etal.[1997[[Lind-[ 
sey & Braun|T997] ). 
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In this review, I present a brief overview of the current 
state-of-the-art in the local helioseismology of magnetic ac- 
tivity, discussing a number of recent findings and some re- 
cent advances in theory that will allow for new insights. 

2 Local Helioseismology of Magnetic Activity 
2.1 Active Region Emergence 

Since its inception, one of the core scientific goals of he- 
lioseismology has been to accurately detect the emergence 
of magnetic activity below the surface. The realization of 
this goal would obviously have considerable implications 
for improving the accuracy of space weather forecasts and 
predictions. There are a number of recent results from lo- 
cal helioseismology that are beginning to show promising 
results on this matter. 

In a series of studies, |Komm et al.| ( |2009a|b| |2007| l and 
|Jain et al.] ( |2012p have used ring diagrams to analyse the sub- 
surface flow systems of a large sample of active and quite- 
Sun regions. Their studies have revealed that flux emer- 
gence is typically accompanied by a number of characteris- 
tics, namely: up-flows, an increase in subsurface zonal flow 
velocity, and an increase in vorticity of subsurface flows, 
while essentially the opposite is observed for decaying flux. 

In another study, |Ilonidis et al. (201 1 1 have used "deep- 
focus" time-distance helioseismology to detect strong 
travel-time reductions 12 - 16 seconds) at great depths 
(~ 42 - 75 Mm), some 1-2 days before high peaks in 
the photospheric magnetic flux rate, for a number of differ- 
ent active regions. These results are somewhat controversial 
however, as the amplitude of these reported travel-time per- 
turbations are significantly larger than theoretical estimates 
based on numerical simulations of emerging flux tubes 
( Birch et al.|^2010) . An independent analysis of the same 



active regions by Braun ( 2012 1, using heUoseismic hologra- 
phy, failed to confirm the reported travel-time anomalies of 
Ilonidis et al. (2011 1, reporting time shifts of order 1 - 2 s 
for the same set of active regions. 

While it has been suggested that subtle differences in 
the measurement geometries/averaging schemes may have 
played a role in the diverging results ( Ilonidis et al.|2012 l, 
the root cause(s) behind why such a large discrepancy ex- 
ists between the two helioseismic methods, and between the 
time-distance results and numerical simulations of emerg- 
ing flux, is currently an open question. Braun ( 20121) makes 
a salient point in that the issue is perhaps best resolved 
through blind tests using artificial data. 

2.2 Sunspot Seismology 

There have been many (at times conflicting) conclusions 
drawn about sunspot and active region structure from local 
helioseismology (see [Gizon et al. |2010| . The predominant 
approach has been to directly interpret helioseismic mea- 
surements and infer the interior properties based on mod- 
els of how perturbations affect these metrics. The physics 



of propagating waves and their interactions with magnetic 
fields are however non-trivial, with mode conversion, which 
is the exchange of energy between the different modes 
of oscillation that are supported by a magnetized plasma 
(fast and slow magnetoacoustic-gravity waves, and Alfven 
waves) appearing to play a key role in helioseismic mea- 
surements. Mode conversion is typically strongest in the re- 
gion where the Alfven wave speed (ca) coincides with the 
sound speed (cj (e.g.|Cally|2005| |2006l|2007 HCrouch et al. 



2005||Schunker & CaUy|2006i r 

The strong magnetic fields which permeate the solar 
surface also produce acute surface phase perturbations that 
obscure helioseismic observations, particularly in sunspot 
penumbrae ( Lindsey & Braun|2003 2005 Schunker et al. 
2005 | l, thereby making magnetic effects much more compli- 
cated and harder to disentangle from other non-linear effects 
such, as radiative transfer effects on Doppler measurements 
( |Rajaguru|20TH|Rajaguru et al.|2007 



Interpreting helioseismic results are often further com- 
plicated when considering the effects of filtering on the data, 
a particularly poignant issue in time-distance helioseismol- 
ogy, where typically either phase-speed filtering or ridge- 
filtering is employed. Phase-speed filtering was originally 
introduced by |Duvall et al.| ( |T997| l in order to improve time- 
distance measurements at small distances. This type of fil- 
tering roughly isolates modes that travel a similar path in 
the solar interior In ridge filtering, modes with a common 
radial order (evident as the "ridges" in a power spectrum) 
are selected ( |DuvaU & Gizon|2000| l. 



In an important study, Braun & Birch (2008i showed 



that the type of filtering employed plays a significant role in 
both the sign and magnitude of measured travel time pertur- 
bations in and around active regions, as phase-speed filters 
can produce anomalous sign-changes as the perturbed ridge 
power falls inside or outside of the filter This in turn has sig- 
nificant implications for the inferred magnitude and lateral 
extent of subsurface structure and flows inversions derived 
from linear inversions of phase-speed filtered helioseismic 
data (see e.g., [Gizon et al.j2009,, Moradi et al.|2010 1. 



2.2.1 Sunspot Flows 

An example of how sensitive inversion results are to the 
type of data filtering employed can be seen in recent in- 
ferences regarding the nature of near-surface flows in and 
around sunspots. In the photosphere, sunspots are typi- 
caUy surrounded by diverging horizontal outflows, termed 
"moat flows". With amplitudes of several hundred m/s, 
they can be seen by local correlation tracking ( |Novem-| 



dBalthasar & Muglach 


2010 


Sobotka & Roudier 2007| Var- 


[gas Dommguez et al. 


2007 


2008 1 and moving magnetic 



features 
[Sheeley 



Brickhouse & Labonte|| 1988) [Gizon et al.||2009t 



969[ ), as well as through Doppler-shift measure- 



ments ( Sheeley 1 1972| l. The moat flow also been detected by 



local helioseismology ( Gizon et al.|2000 1. 
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In a recent result, the horizontal and vertical extent of 
the near-surface flows associated with the sunspot in AR 
9787 was characterised by |Gizon et al.|p009| l. Using lin- 
ear inversions of ridge-filtered time-distance travel times 
(in conjunction with Born-approximation sensitivity kernels 



of Birch & Gizon (2007 1), they detected a strong outflow 
of several hundred m/s extending beyond the penumbra of 
the sunspot, down to a depth of around 4.5 Mm. These 
results were later confirmed by Moradi et al. ( 2010| l us- 
ing high resolution ring diagram analysis. More recent 3D 
ring-diagram inversions by Featherstone et al. ( 201 l|l (u s- 
ing the Born-approximation kernels of Birch et al. ( 2007 1), 
and inversions using helioseismic holography ( Braun et al.| 
|201 1[ ) have also confirmed similar near-surface outflows as- 
sociated with sunspots. These results contrast strongly with 
those of [Zhao et al.| ( [200 1'), who used linear inversions 
of phase-speed filtered time-distance travel times and ray- 
approximation sensitivity kernels, to infer near-surface in- 
flows (i.e., flows in the opposite direction to the observed 
moat flow) associated with sunspots. 



2.2.2 Sunspot Structure 

Almost all of the various local helioseismic diagnostic tools 
have also been used to infer changes in the wave speed 



in sunspots (e.g., Baldner et al. 



ICouvidat et"aL|[2006,; Fan et al. 



2011 



1995 



Basu et al. 2004 



Kosovichev et al.| 



of work on this topic, there is still no general agreement 
on the subsurface structure of sunspots. Gizon et al. (2009) 
and [Moradi et al.| ( |2010[ ) have shown that inferences made 
by time-distance and ring-diagram analysis with regards to 
subsurface wave-speed structure are surprisingly contradic- 
tory, with the inversions producing subsurface wave-speed 
profiles of opposite signs and different amplitudes. 

The fact that neither method fully accounts for the de- 
tails of the measurement procedures (e.g., the effects of 
phase-speed filtering), and that both inversion methods use 
sensitivity functions that do not explicitly include the direct 
effects of the magnetic field, are perhaps the most poignant 
factors which contribute to such stark disagreements be- 
tween the two methods. However, when ,Moradi et al.| ( |2010| l 
and Gizon et al. ( 2010| l compared the linear inversion re- 
sults with a number of phenomenological and purely nu- 
merical models of subsurface wave-speed structure (e.g., 
[Cameron et aL][20ni [Crouch et al.||2005l [Fan et al.|[l^ 



Rempel et al. [[2009b i, they found that the ring-diagram in- 



versions and the phenomenological/numerical models were 
actually in agreement - to an extent - displaying positive 
wave-speed perturbations down to a depth of ~ 2 Mm. 

As solar imaging hardware becomes increasingly so- 
phisticated, with the Helioseismic and Magnetic Imager 



(HMI; Scherrer et al.|2 012 ) onboard SDO providing contin- 
uous Doppler and magnetic images covering the full solar 
disk with fourfold better spatial resolution than its prede- 
cessor instrument (SoHO-MDI), the need for refined numer- 



ical modelling/computational techniques becomes apparent 
to realize the full potential that these observations offer 

3 Numerical Modelling 
3.1 Forward Modelling 

In the last few years, several advanced 3D (mainly linear) 
parallel codes have been developed to conduct computa- 
tional simulations of magnetohydrodynamic (MHD) waves 
in sunspots, where strong magnetic fields and severe grav- 
itational stratification produce complex behaviours in the 
waves we use to probe the interior (e.g., [Cameron et al. 



2008; Hanasoge 2008; Khomenko et al."2009| [Parchevsky 
&' Kosovichev 2009, Shelyag et al. 2009). Through the 



process known as forward modelling - constructing com- 
putational models that mimic wave propagation through 
sunspots and matching the resultant wave statistics with ob- 
servations - we can refine the helioseismic inferences that 
are currently fooled by these wave interaction and conver- 
sion processes. The major advantage of this approach is that 
it permits the study of wave propagation by providing the 
freedom to choose various background models, without the 
complications of solving for convection. This substantially 
reduces the computational expense for detailed statistical 
studies. Interpretations of available results from the SoHO 
and SDO missions can thus be re-evaluated with these stud- 



2000 Zharkov et al. 2007 1. But despite the long history ies. 



A successful example of this approach is from the anal- 
ysis undertaken by Cameron et al. ( 201 l| l, who model the 
propagation of /, pi and p2 planar wave packets through 
a shallow, semi-empirical, non-magnetohydrostatic (MHS) 
sunspot model of AR 9787. In order to assess the valid- 
ity of the underlying sunspot model, the authors compared 
the simulated wave field to the observed SoHO-MDI cross- 
covariance functions. They found a good quantitative match 
between the phase shifts and amplitudes of the simulations 
and actual observations of AR 9787, which is strong evi- 
dence that the wave signal is dominated by the near-surface 
layer of the sunspot. 

A recent (observational) study of sunspot-induced 
acoustic wave scattering by Zhao et aL] ( [2011 1 also tends to 
favour this assertion. Using data from the HMI instrument, 
they analysed the wavefunctions of scattered waves with the 
incident waves for the / to mode for two sunspots. They 
observed that the magnitude of the scattered wave relative 
to that of the incident wave, decreases from the /- to the 
P5 mode, indicating that the region producing the scattered 
wave is shallower than the depth of the /-modes. 

3.2 Inverse Modelling 

The next step in the numerical modelling process is to iden- 
tify a clear and efficient way to alter sunspot models and the 
flows around them in order to better match observations. 
This will entail the solving of an inverse problem, which for 
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local helioseismology is to use measured helioseismic sig- 
natures (e.g., travel times) to estimate physical conditions 
in the solar interior. While recent progress has been made 
in improving the tradition al linear sensitivity kernels and 
inversion algorithms (e.g., Jackiewicz et al. 2012 Svanda 



et al. 201 1 1, these methods have traditionally relied on back- 



ground models being translationally invariant (e.g. Gizon 



|& Birch||2002[ ), which is problematic (as discussed previ- 
ously), since the nature of waves in solar active regions 
is significantly altered as they are converted into magneto- 
acoustic waves. Given the mathematical difficulties in mod- 
elUng such complex phenomena, current linear inversions 
do not take the effects of the sunspot magnetic field into ac- 
count, nor the possibility that such wide deviations from the 
quiet Sun introduce a non-linear relation between observed 
travel-time shifts and causative agents. 

A significant breakthrough has been made on this 
front however, with the development of an improved in- 
verse method for helioseismology. The adjoint method. 



employed by Hanasoge et al. (2012 20111, is a form of 
partial-differential-equation constrained optimization that 
describes a means of extracting derivatives of a misfit func- 
tion with respect to model parameters through finite compu- 
tation. The technique allows for the computation of sensi- 
tivity kernels around arbitrary background models, thereby 
paving the way for non-linear iterative inversions. On an- 
other front, [CFouch et aL] ( |2011| l have managed to success- 
fully invert for the sound-speed and magnetic field in a sim- 
ple ID model. 

3.3 Numerical Issues 

While advances in numerical modelling techniques has al- 
lowed us to make some significant inroads into a number of 
outstanding issues in the local helioseismology of magnetic 
activity, they are not without their own problems. In this 
section, I discuss a number of the more pressing numerical 
issues, affecting both forward and inverse modelling, that 
need to be addressed in order for us to have more confidence 
in the numerical simulations. 

i) Convective instability: As the Sun is highly convec- 
tively unstable in the outer 2-3% layers, simulating small 
perturbations around this state is not feasible (a particu- 
larly pertinent issue for the linear wave propagation codes). 
A somewhat passive approach which has been commonly 
adopted by the helioseismology community, is to slightly al- 



ter the background model (typically Model S; Christensen 



convective instability associated with the outermost layers 
currently has not been explored. 

ii) Sunspot models: Numerical models of sunspots are 
a necessity for both forward and inverse modelling. As de- 
( 2010| l, there is certainly no short- 



tailed in Moradi et al. 



age of models available for helioseismologists to play with, 
which in itself is not a disadvantage, as examining a num- 
ber of different models is essential to test the robustness of 
helioseismic inferences. Ideally, sunspot models should be 
constrained by surface observations and constructed with 
the subphotospheric layers in force balance and the atmo- 
spheric region force-free. 

On the other hand, fully realistic magnetoconvective 
simulations of sunspot structure are now a reaUty (e.g.. 



Heinemann et al. 2007 Rempel 2012 Rempel et al. 



2009a|b I, but currently too cost prohibitive to be employed 



for the exhaustive statistical studies which are inherently 
required for computational helioseismology. Nonetheless, 
such simulations are a valuable tool for improving both for- 
ward and inverse modelling techniques. A recent analysis 



by Braun et al. ( 20 12) has shown a very good quantitative 
agreement between the travel-time shifts of a magnetocon- 
vective simulation and those of two actual sunspots mea- 
sured using helioseismic holography. 

iii) Alfven wave speed: One of the biggest difficulties in 
simulating 3D MHD wave propagation in the Sun is due to 
the excessively large value that Ca attains above the surface, 
approaching values of several hundred km/s. This causes the 
wavelengths of both the fast and Alfven waves to become 
rather large, resulting in an extremely stiff numerical prob- 
lem for the simulation codes, with the simulation timestep 
(At) being highly constrained by the Courant-Friedrichs- 
Lewy (CFL) condition (typically A< ^ Az/ca; where Az 
denotes the vertical grid resolution). As such, in order to 
simulate artificial data sets on the time scales required for 
computational helioseismology (in a finite amount of time), 
a number of wave propagation codes apply somewhat sim- 
plistic approximations to treat the overlying model-sunspot 
atmosphere (e.g., [Braun et aL 2012 Cameron et al.||20lT] 



Hanasoge]|2008 Rempel et al. 



2009b I. The most common 



approach has been to apply a "limiter" to moderate the 
action of the Lorentz forces when Csica becomes unreason- 
ably big, with Ca being typically (smoothly) capped at val- 
ues of between 20-60 km/s above the surface. 

The physical implications of artificially limiting on 
the seismology are not well understood however. While 
some choices for these Umiters have been discussed be- 



[Dalsgaard et al.| 1996[) in order to keep the Brunt- Vasala fre- fore (e.g., [Braun et al.||2012| [Cameron et aL[[2008l [20TT 



quency positive ( Hanasoge et al.|2006 Parchevsky & Koso 



vichev|2007[[Schunker et al.|2011[[Shelyag et al.|2009[ ). 

While this method produces solar-like power spectrums, 
it is far from an ideal solution, as changing the background 
in such a manner effectively shifts the eigenfunctions and 
eigenfrequencies of the resonant modes, and may also af- 
fect the seismic reciprocity of the system ( [Dahlen & Tromp| 
1998 1. A clear idea of how to successfully deal with the 



[Rempel et al . 2009b), no serious effort has yet gone into 
understanding their impact on the near- and far-scattered 
wave-field. This is a particularly pertinent issue, given the 
critical role that the exponentially increasing Ca appears to 
play in the fast-Alfven wave mode conversion process in 
sunspot atmospheres (Cally|2012[[Cally & Goossens|2008| l. 
In a series of theoretical studies using the cold plasma ap- 
proximation ( [Cally & Andries|2010|[Cally & Hansen|201 1 
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Hansen & Cally|2012 Hanson & Cally'201 and more re 
alistic model sunspot atmospheres ( Felipe,2012 Khomenko 
|& Cally||2012| i, it has been demonstrated that conversion 
to Alfven waves occurs at, and beyond, the fast wave re- 
flection height (located approximately where Ca ^ oj/kh', 
where lu denotes the wave frequency and kh the horizon- 
tal wavenumber), being spread over many scale heights for 
wavenumbers typical of local helioseismology. The fast- 
Alfven conversion process appears to be most efficient for 9 
(field inclinations from vertical) ^ 30° - 40°, and (f) (the 
angle between the magnetic field and wave propagation 
planes) -60° -80°. 

So this leads to a very important question: What are the 
implications of the returning fast and Alfven waves for the 
seismology of the photosphere? While the answer to this 
question currently remains elusive, it does appear that lim- 
iting Ca, in the manner that we have been doing, could be 
potentially problematic for active region seismology. 

4 Conclusion 

It is perhaps not too presumptuous to say that local helio- 
seismology is the greatest benefactor from improved obser- 
vations provided by SDO. Since 2010, the HMI instrument 
not only provides higher-resolution Doppler and continuum 
images of the solar disk, but also high quality line-of-sight 
and vector magnetograms - the latter are extremely valuable 
resources for helioseismologists. Reliable measurements of 
the photospheric magnetic field are the key to constrain- 
ing the near-surface layers of MHS sunspot models, as de- 
tailed models of the surface layers are a necessity in order to 
probe the deeper structure of sunspots. Whether the "mono- 
Hthic " (e.g., |Cowling]|1946| |1957| |1976|), or "spa ghetti" 
(e.g., |Parker| 19751 1 1 979| |Spruit| 198 1 [ |Zwaan| 1 98 1 | l model 
best explains their subsurface magnetic field structure, still 
remains to be seen. 

But as we wait for new observational results to roll in, 
we should, in the meantime, continue to improve the ac- 
curacy and efficiency of our numerical forward modelling 
techniques to address the inconsistencies in the helioseis- 
mic inferences of both flux emergence and active region 
seismology studies. Both areas are also in pressing need of 
detailed parametric studies, the results from which can be 
utiUsed to fine-tune translationally-variant sensitivity ker- 



nels ( Hanasoge et al.||2012) , required for any future non- 
hnear inversions of subsurface sunspot structure and dy- 



namics. 
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